Orbital Kondo effect in Cobalt-Benzene sandwich molecules 
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We study a Co-benzene sandwich molecule bridging the tips of a Cu nanocontact as a realistic 
model of correlated molecular transport. To this end we employ a recently developed method for 
calculating the correlated electronic structure and transport properties of nanoscopic conductors. 
When the molecule is slightly compressed by the tips of the nanocontact the dynamic correlations 
originating from the strongly interacting Co 3d shell give rise to an orbital Kondo effect while the 
usual spin Kondo effect is suppressed due to Hund's rule coupling. This non-trivial Kondo effect 
produces a sharp and temperature-dependent Abrikosov-Suhl resonance in the spectral function at 
the Fermi level and a corresponding Fano line shape in the low bias conductance. 



The discovery of ferrocene and bisbenzene chromium 
[l[ over half a century ago was the starting point 
for experimental and theoretical work both in chem- 
istry and physics concerning the intricate properties of 
organometallic compounds. The investigation of these 
and related sandwich complexes is driven by their rel- 
evance in various chemical applications (e.g. cataly- 
sis) and more recently also because of their prospective 
nanotechnological applications for example as molecular 
magnets Q or spintronic devices 0. 

Molecules with a transition metal (TM) center coupled 
to aromatic groups are also of high interest from a funda- 
mental point of view. The strong electronic correlations 
in the d shell of the TM can modify the ground state and 
electronic transport properties of such molecules, lead- 
ing to many-body phenomena like the Kondo effect |4[ 
as recently observed in TM-phthalocyanine molecules [5| . 
The importance of dynamical correlations in nanoscopic 
devices in general is further substantiated by the recent 
observation of the Kondo effect in nanocontacts made 
from TMs @, 0] • Also recently a so-called underscreened 
Kondo effect has been reported for a molecule trapped 
in a breakjunction Q. While the Kondo effect is com- 
monly associated with the screening of a local magnetic 
moment by the conduction electrons, it is also possible 
that another internal degree of freedom associated with a 
degeneracy gives rise to a Kondo effect 0] . One example 
is the so-called orbital Kondo effect where the pseudo- 
spin arising from an orbital degeneracy is screened by 
the conduction electrons ■ 

Such complex many-body phenomena call for a theo- 
retical description beyond the standard treatment with 
Kohn-Sham (KS) density functional theory (DFT) which 
cannot capture many-body physics beyond the effective 
mean-field picture. The incorrect behavior of the KS 
DFT for strongly correlated systems can be remedied by 
augmenting the DFT with a local Hubbard-like interac- 
tion. This DFT++ approach is the de facto standard in 
the theory of solids |11| . Recently this approach has been 
adapted to the case of nanoscopic conductors @, [H, EH ■ 

In this letter we apply this method to investigate 



the transport properties of cobalt benzene sandwich 
molecules (CgH^Co^eHe) (C0BZ2 in the following). We 
find that dynamical correlations in the Co 3c? shell give 
rise to an orbital Kondo effect in a doubly degenerate 
level of the Co 3d shell while the usual spin Kondo effect 
is suppressed due to Hund's rule coupling. 

We perform DFT calculations using the CRYSTAL06 
code Q employing the LDA 0, PW91 [HI and the 
hybrid functional B3LYP [13], and using the all electron 
Gaussian 6-31G basis set. The geometries of the wires 
were relaxed beforehand and kept fixed during the cal- 
culations. The geometry of the molecule in contact with 
the wires was relaxed employing the B3LYP functional 

In order to capture many-body effects beyond the 
DFT level, we have applied the DFT+OCA (Density 
Functional Theory + One-Crossing Approximation) for 
nanoscopic conductors developed by one of us in earlier 
work @,im • To this end the system is divided into three 
parts: The two semi-infinite leads L and R and the de- 
vice region D containing the molecule and part of the 
leads. The KS Green's function (GF) of region D can 
now be obtained from the DFT electronic structure as 
G° (w) = (w + n - H% - E L (w) - ErM)- 1 where H° 
is the KS Hamiltonian of region D and Sl ) r(w) are the 
so-called lead self-energies which describe the coupling 
of D to L and R and which are obtained from the DFT 
electronic structure of the nanowire leads. 

Next the mean-field KS Hamiltonian is augmented by 
a Hubbard-like interaction term J2ijki Uijud\djdidk that 
accounts for the strongly interacting electrons of the Co 
3d shell. Here we use a simplified interaction which 
only takes into account the direct Coulomb repulsion 
U = Uijij and the Hund's rule coupling J = Uijji. These 
are different from the bare interactions due to screen- 
ing processes. Here we assume that the interactions are 
somewhat increased as compared to their bulk values 
since the screening should be weaker than in bulk. Hence 
we use U = 5 eV and J = 1 eV; indeed we find that our 
results are qualitatively stable for a reasonable range of 
values for U and J (see below). 
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The interacting Co 3d shell coupled to the rest of the 
system (benzene+leads) constitutes a so-called Ander- 
son impurity model (AIM). The AIM is completely de- 
fined by the interaction parameters U and J, the en- 
ergy levels td of the 3d orbitals and the so-called hy- 
bridization function A^(a;). The latter describes the 
(dynamic) coupling of the Co 3d-shell to the rest of 
the system and can be obtained from the KS GF as 
Ad(w) = ui + jU — e d — [G"(cj)] _1 where fj, is the chemical 
potential, e° are the KS energy levels of the 3d orbitals 
and G d (uj) is the KS GF projected onto the 3d subspace. 
The energy levels €d are obtained from the static crystal 
field e d where as usual in DFT++ approaches a double 
counting correction (DCC) has to be subtracted to com- 
pensate for the overcounting of interaction terms. Here 
we employ the so-called fully localized or atomic limit 
DCC [19\E dc = U- (N 3d -±)-J (^f -J). 

The AIM is solved within the OCA [ill HI- This 
yields the electronic self-energy Ed(^) which accounts 
for the electronic correlations of the 3ci-electrons due 
to strong electron-electron interactions. The correlated 
3d GF is then given by G d = ([G°]- x - E d + E dc )-\ 
Correspondingly, the correlated GF for D is given by 



G D = ([G D _ 



Edc) where and Ed c only 



act within the 3d subspace. From Gd we can calculate 
the transmission function T(u>) = Ti'IFlGqFjjGd] where 
T a = i(S Q — EjJ. For small bias voltages V, the trans- 
mission yields the conductance: Q(V) = (2e 2 / h)T (eV) . 

We consider a single C0BZ2 molecule in contact with 
two semi-infinite Cu nanowires as shown in Fig. QJi. The 
C0BZ2 molecule is the smallest instance of a general 
class of M„Bz n _|_i complexes, where M stands for a metal 
atom, that have been prepared and investigated 0, 21 



The semi-infinite Cu wires exhibit the hexagonal sym- 
metry of the molecule and correspond to the (6,0) wires 
described in Ref. [22J. We investigate the system at three 
different Cu-tip-Co distances 3.6 A, 4.0 A and 4.3 A (see 
Fig. [lji,b). As can be seen from Fig. [1J> the Bz-Bz dis- 
tance h varies between 3.4 A and 3.65 A depending on 
the distance d of the Cu tip to the Co atom in the center 
of the molecule. At distances d — 3.6 A and d = 4.0 A 
the molecule is slightly compressed compared to its free 
height of about h = 3.6 A, whereas it is slightly stretched 
at d = 4.3 A. 

The hexagonal symmetry of the system leads to a lift- 
ing of the degeneracy of the five 3d orbitals. The sym- 
metry adapted representations are: the A\ group con- 
sisting of the d^ z 2_ r 2 orbital only, the doubly degenerate 
Ei group consisting of the d xz and d yz orbitals and the 
E2 group consisting of the d xy and d x i_ y i orbitals. 

Figs. QJ+d show the imaginary parts of the hybridiza- 
tion functions A<j(w) calculated from the LDA electronic 
structure. The imaginary part of the hybridization func- 
tion exhibits a distinct peak close to the Fermi level (Ep) 
in the E2 channel, whose position, width and height de- 
pend significantly on the molecular geometry, specifically 
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Figure 1. (Color online) (a) Atomic structure of C0BZ2 
molecule in a Cu nanocontact. (b) Size h of C0BZ2 molecule 
as a function of the distance d between Co atom and Cu tip 
atoms, (c) Orbitally resolved imaginary part of the hybridiza- 
tion function for d = 3.6A. (d) Total hybridization function 
(all Co 3d orbitals) for different d. 



on the Bz-Co distance. The other channels E\ and A\ 
show only a negligible hybridization close to Ep. The 
dominant feature stems, similarly as shown for graphene 
12 , [23| from hybridization with the 7r 2 orbital state of 
the benzene rings. The feature does not depend quali- 
tatively on the DFT functional used, as we have found 
the same feature within GGA and also in B3LYP calcu- 
lations. The presence of strong molecular resonances in 
the hybridization function makes this case different from 
the case of nanocontacts with magnetic impurities where 
the hybridization functions arc generally much smoother 
(see Ref. Q for comparison). 

Fig. [2k compares the correlated spectral functions 
Ad(uj) ~ ImGd of the Co 3d electrons for the three 
distances considered here at high temperatures on a 
large energy scale. The spectra vary considerably as 
the distance d changes. Most importantly, for d around 
3.6 A when the molecule is slightly compressed, a sharp 
temperature-dependent peak appears right at Ep, as can 
be seen from Fig. [2Jd. The peak is strongly renormalized 
(i.e. it only carries a small fraction of the spectral weight) 
due to the strong electron-electron interactions. 

The sharp peak in the spectral function at Ep that 
starts to develop already at temperatures of kT 
f J .01 eVw 120 K stems from the E2 channel which is the 
only channel with appreciable hybridization near Ep, see 
Figs. [TJd+c. Correspondingly, the transmission function 
(Fig. [2}:) shows a Fano-like feature around zero energy. 
A renormalized, sharp and temperature-dependent res- 
onance in the spectral function at Ep is commonly as- 
sociated with the Kondo effect. Looking at the orbital 
occupations, we find that the E2 channel that gives rise 
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Figure 2. (Color online) (a) Spectral functions of Co 3d shell 
for different d at temperature T ~1200 K. (b,c) Spectral and 
transmission functions of the molecule at d = 3.6 A for differ- 
ent temperatures. 



to the resonance for d = 3.6 A has an occupation of 
about 2.8 while the total occupation of the 3d shell is 
Nd ~ 7.5. The fractional occupation numbers indicate 
the presence of valence fluctuations where the charges in 
the individual impurity levels fluctuate in contrast to the 
pure Kondo regime where these fluctuations are frozen. 

Analyzing the atomic states of the Co 3d-shell con- 
tributing to the ground state of the system we find that 
the principal contribution (~ 45%) is an atomic state 
with 8 electrons and a total spin of S = 1 (d 8 , S = 1) as 
shown in Fig. [3Ji. The total spin 1 stems from holes in 
the E2 and A\ channels. The charge fluctuations in the 
Ei channel are mainly due to the contribution (~ 17%) 
of an atomic (d 7 , S = 3/2) state. There are considerably 
weaker contributions (~ 4%) from atomic (d 7 ,S = 1/2) 
and (d 9 , S — 1/2) states. The individual contributions of 
the remaining atomic states are very small (below 1%) 
but add up to a total contribution of 34%. 

By exclusion of individual atomic states from the cal- 
culation of the spectra we can determine which fluctu- 
ations are responsible for the different spectral features. 
We find that the fluctuations between the (d 8 ,S = 1) 
and the (d 7 ,S = 3/2) states are primarily responsi- 
ble for the three spectral features close to Ep including 
the sharp Kondo-like peak right at Ef, as illustrated in 
Fig- Eh- Also the broad peak around 4 eV below Ep orig- 
inates from these fluctuations while the broad peak about 
2.5 eV above Ep arises from fluctuations between the 
(d 9 ,S = 1/2) atomic state and the principal d 8 atomic 
state. The two peaks in the spectral function nearest to 
the Kondo resonance arise from the strong energy depen- 
dence of the hybridization function whose real part has 
poles just below and above Ep roughly at the positions 
of these two spectral features. 



Note that the fluctuations from the (d 8 ,S = 1) to the 
(d 7 ,S = 3/2) states that give rise to the Kondo-like peak 
at Ep actually cannot lead to a spin Kondo effect since 
the spin 3/2 of the d 7 state is higher than the spin 1 of 
the principal d 8 state. Instead the fluctuations between 
the (d 8 ,S = 1) and the (d 7 ,S = 3/2) states which give 
rise to the Kondo-like resonance at Ep correspond to an 
orbital Kondo effect in the doubly-degenerate E2 levels 
of the Co 3d shell as illustrated in the upper panel of 
Fig. [3t. Here the index labeling the two orbitals with 
Ei symmetry takes over the role of a pseudo spin. In 
the principal d 8 atomic state the E2 levels are occupied 
with three electrons and hence have a pseudo spin of 
1/2. By excitation to the (d 7 , S = 3/2) state the electron 
with minority real spin and with some pseudo spin state 
is annihilated. By relaxation to the principal electronic 
(d 8 ,S =1) state a minority real spin electron can now 
be created in one of the two pseudo spin states. Those 
processes that lead to a flip of the pseudo spin then give 
rise to the orbital Kondo effect and the formation of the 
Kondo peak at Ep. 

The absence of a normal spin Kondo effect where the 
total spin 1 of the principal d 8 atomic state is screened, 
is easily understood on the following grounds: First, in 
general the Kondo scale decreases exponentially with in- 
creasing spin of the magnetic impurity [24( . In addition, 
here the A\ level does not couple at all to the conduction 
electrons around Ep (no hybridization). Thus the spin 
1/2 associated with it cannot be flipped directly through 
hopping processes with the conduction electron bath. 

On the other hand, an underscreened Kondo effect as 
reported in Ref. where only the spin 1/2 within the Ei 
shell is screened is also suppressed compared to the or- 
bital Kondo effect due to Hund's rule coupling: Screening 
of the spin 1/2 in the E2 shell can take place by fluctu- 
ations to the (d 7 , S — 1/2) state. However, the Hund's 
rule coupling J favors the high spin (d 7 , S = 3/2) state 
over the low spin (d 7 , S = 1/2) state as can also be seen 
from the smaller weight of the latter compared to the for- 
mer. Hence the Kondo scale is considerably lower for the 
underscreened Kondo effect than for the orbital Kondo 
effect found here. At lower temperatures (not accessible 
with OCA) the two Kondo effects may in fact coexist. 

We have also checked the dependence of the 
LDA+OCA spectra on the DCC as well as on the inter- 
action parameters U and J (not shown [IBj]). In general 
we find the spectra and also the Kondo peak to be qual- 
itatively robust against shifts of the impurity levels in 
energy over a range of several electron volts, and changes 
of U between 3 and 7 eV, and of J between 0.6 to 1 eV. 
As expected for the Kondo effect the sharp resonance 
stays pinned to the Fermi level when shifting the impu- 
rity levels in energy, and only height and width somewhat 
change. 

Stretching the molecule by displacing the tips of the 
Cu nanowires the Kondo resonance and the concomitant 
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Figure 3. (Color online) (a) Fluctuations between the atomic 
states that give rise to the spectral features shown in (b). (b) 
Spectral function around Ef for d = 3.6 A. Arrows indicate 
spectral features arising from the fluctuations between atomic 
states shown in (a), (c) Orbital (red) and spin flip (blue) fluc- 
tuations in the E2 channel that can lead to an orbital Kondo 
effect and to an underscreened spin Kondo effect, respectively. 



Fano line shape in the transmission disappear for dis- 
tances d > 4 A (not shown). This is accompanied by 
an increase of the occupation of the Co 3d shell. The 
new regime is characterized by a strong valence mixing 
between the d 8 and d 9 atomic state of roughly equal con- 
tribution indicating that the system is now in the so- 
called mixed valence regime (see e.g. Ref. 0, Chap. 5). 
Hence the orbital Kondo effect and the associated spec- 
tral features can be controlled by stretching or compress- 
ing the molecule via the tip atoms of the Cu nanocon- 
tact. This strong sensitivity on the molecular conforma- 
tion stems from the sharp features in the hybridization 
function which change considerably when the molecule 
is stretched or compressed as can be seen from Fig. [TU. 
This peculiar behavior is qualitatively different from the 
case of the nanocontacts containing magnetic impurities 
where the hybridization functions are much smoother 0] . 

In conclusion, we have shown for a C0BZ2 sandwich 
molecule coupled to Cu nanowires that the dynamic cor- 
relations originating from the strongly interacting Co 3<i 
electrons give rise to an orbital Kondo effect in the dou- 
bly degenerate E2 levels while a spin Kondo effect is sup- 
pressed by Hund's rule coupling. Due to the sensitivity 
of the electronic correlations on the molecular confor- 
mation it is possible to control the appearance of the 
orbital Kondo effect by stretching and compressing the 
molecule with the tips of the Cu leads. It should be pos- 
sible to prepare a setup similar to the one considered here 
in an actual experiment and measure the orbital Kondo 
effect, e.g. by contacting a CoBz 2 molecule deposited on 
a Cu(lll) surface with a scanning tunneling microscope. 
Finally, we mention that we have also explored the case 



of a NiBz2 sandwich. In this case we do not find a Kondo 
effect. Instead the system is in the mixed valence regime 
characterized by a strong peak with weak temperature 
dependence close to but not at Ep [25| . 
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